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Abstract. We have derived energy conservation equations from the quaternionic 
Newton's law that is compatible with Lorentz transformation. This Newton's law yields 
directly the Euler equation and other equations governing the fluid motion. With this 
formalism, the pressure contributes positively to the dynamics of the system in the 
same way mass does. Hydrodynamic equations are derived from Maxwell's equations 
by adopting an electromagnetohydrodynamics (EMH) analogy. In this analogy the 
hydroelectric field is related to the local acceleration of the fluid and the Lorentz 
gauge is related to the incompressible fluid condition. An analogous Lorentz gauge in 
hydrodynamics is proposed. We have shown that the vorticity of the fluid is developed 
whenever the particle local acceleration of the fluid deviates from the velocity direction. 
We have also shown that Lorentz force in electromagnetism corresponds to Euler force 
in fluids. Moreover, we have obtained Gauss's, Faraday's and Ampere's -like laws in 
Hydrodynamics. 



PACS numbers: 41.20. Jb; 47.10.-g; 47.10.A- 
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1. Introduction 

We have recently formulated Maxwell equations using quaternion (Arbab and Satti, 
2009 & Arbab, 2009). In this formalism, we have shown that Maxwell equations have 
not included the magnetic field produced by the charged particle. Consequently, we have 
also shown that this magnetic field is given by an equation equivalent to Biot-Savart 
law. Moreover, we have shown that the magnetic field created by the charged particle is 
always perpendicular to the particle direction of motion. We have then introduced the 
quaternionic continuity equation and compared it with the ordinary continuity equation. 
Consequently, we have found that Maxwell equations predict that the electromagnetic 
fields propagate in vacuum or a charged medium with the same speed of light in vacuum 
if this quaternioic continuity equation is valid. 

In the present paper we explore the idea of formulating the fluid equations using 
the quaternion description. Owing to the previous analysis done in electromagnetism, 
we would like to further establish an analogy between quaternionic hydrodynamic and 
quaternionic Maxwell's equations. Accordingly, it turns out that the quaternionic 
Newton's equation gives directly the non-relativistic limit of the energy momentum 
conservation equation. These equations are the Euler equation, the continuity equation, 
the energy momentum conservation equation. Thus, an analogy between hydrodynamics 
and electrodynamics is shown to exist. This analogy guarantees that one can derive the 
appropriate equations of the former system from the latter or vice versa. This analogy 
is quite impressing since it allows us to visualize the flow of the electromagnetic field as 
resembling the flow of a fluid. A quaternionic Newton's second law is developed and a 
generalized Newton's law of motion is then obtained. The application of the generalized 
continuity equations (GCE) to the generalized Newton's law yields the non-relativistic 
equations governing the motion of fluids. It is interesting to notice that generalized 
continuity equations are Lorentz invariant. The fluid vorticity is found to arise whenever 
the local acceleration of the fluid particles deviates from the velocity direction. The 
physical properties and the equations governing the hydrodynamics are derived. We 
remark here that these equations can be put in in analogous forms to Gauss's, Faraday's 
and Ampere's laws. Moreover, we have found that the diffusion equation is compatible 
with the GCEs. We therefore, emphasize that the GCEs should append any model 
dealing with fluid motion. According to this analogy and since the electrodynamics is 
written in terms of the electric and magnetic fields, E e and B respectively, we should 
write the hydrodynamics equations in terms of the hydroelectric field and vorticity 
Cu and will then obtained analogous equations to Maxwell equations. This is legitimate 
because of the analogy that exists between the two paradigms permits, viz., Eh E e 
and u B, v A and x = y ^ where E h = — 1|, u = V x v . 
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2. The Continuity equation 

We have recently explored the application of quaternions to Maxwell's equations (Arbab 
and Satti, 2009). We have found that the quaternionic Maxwell's equation reduces to 
the ordinary Maxwell's equations but predicts the existence of a scalar wave competing 
with the electromagnetic field traveling at the speed of light in vacuum. Quaternions are 
found to have interesting properties. The multiplication rule for the two quaternions, 
A = (oq , A) and B = (b , B) is given by (Tait, 1904) 

AB = (a b - A - B , a B + Ab + A X B) . (1) 

Therefore, the ordinary continuity equation is transformed into a quaternionic continuity 
equation as 



VJ 



0, (2) 



where 



i|,v), 7=(ipc,f). (3) 
The scalar and vector parts of Eq.(2) imply that 

V-/ + |=0, (4) 

Vp + if = 0, (5) 

and 

VxJ = 0. (6) 

We call Eqs.(4)-(6) the generalized continuity equations (GCEs). In a covariant form, 
Eqs.(5)-(6) read 

V = 0, d„J v -d v J ll = Q. (7) 

Hence the GCEs are Lorentz invariant. We remark that the GCEs are applicable to any 
flow whether created by charged particles or neutral ones. 

3. Newton's second law of motion: 

The motion of the mass (m) is governed by the Newton's second law. The quaternionic 
Newton force reads (Arbab and Satti, 2009) 

F = -mV{VV) , (8) 

where 

Z -P,F), V={ic,v). (9) 
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The vector part of Eq. (8) yields the two equations 

dv 

and 



™(^ + V ( V — ) -v x (V x v) ) , 



-l _ v dv 

Vx„ = --x-. (11) 

The scalar part of Eq.(8) yields the two equations 

v dv 

and 

tf. x i?) = 0. (13) 
For a continuous medium (fluid) containing a volume 1/, one can write Eq.(10) as 



P = mc U V-v + --^ \ , 



,v 
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p(^ + VC-)-v*x(Vxv))=f, (14) 

where m = p V and / = y . Using the vector identity |V(u -v) = iJx(Vxi;) + (fV)i), 
Eq.(14) becomes 

p(j% + (v-V)v > j = /. (15) 

This is the familiar Euler equation describing the motion of a fluid. For a fluid moving 
under pressure (P r ), one can write the pressure force density as 

f P = -VP r , (16) 

so that Eq.(15) becomes 

(dv - \ 

Using Eq.(5), Eq.(12) can be written as 

nil —>—>—> — # 

_ + V .5 = /.ir j 5=(pc 2 ) ? 7, M = pt; 2 , P = /.iT. (18) 

Eq.(18) is an energy conservation equation, where S is the energy flux and u is the 
energy density of the moving fluid. With pressure term only, Eq.(18) yields 
dv -* -* 

— + V-(pc 2 + P r )v = P r V-v. (19) 

The source term on the right hand side in the above equation is related to the work 
needed to expand the fluid. It is shown by Lima et al. (1997) that such a term has 
to be added to the usual equation of fluid dynamics to account for the work related 
to the local expansion of the fluid. It is thus remarkable we derive the fundamental 
hydrodynamics equations from just two simple quaternionic equations, the continuity 
and Newton's equation. For incompressible fluids V • v = so that Eq.(19) becomes 

^ + V-(pc 2 + P r )v = 0, (20) 
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which states that the pressure contributes positively to the energy flow. This means the 
total energy flow of the moving fluid is 

Stotai = (pc 2 + P r ) v, (21) 

and the total momentum density of the flow is given by 

p = (p+^>7, (22) 
cr 

which must be conserved. This is analogous to the general theory of relativity where 
the pressure and mass are sources of gravitation. Equation (20) states also that there 
is no loss or gain of energy for incompressible fluids. However, when viscous terms 
considered loss of energy into friction will arise. In standard cosmology the general 
trend of introducing the bulk viscosity (r/j*| is by replacing the pressure term P r by the 
effective pressure (Weinberg, 1972) 

Pes. =P r -vV-v. (23) 

Substituting this in Eq.(17) and defining the vorticity of the fluid by u = V x v, we get 



fdv - \ 

p ( — + (v ■ V) v j = -VP r + r/V\7 + r/V x cD , (16a) 

which reduces to the Navier-Stokes equation for irrotational flow (u = 0). Equation 
(20) can be put in a covariant form as 

d^ v = , T fiu = (p+^) v^v u - P r g, w . (24) 

where T^ u is the energy momentum tensor of a perfect fluid, u M is its velocity and is 
the metric tensor with signature (+++-). It is interesting to remark that we pass from 
quaternion Newton's law to relativity without any offsetting. This is unlike the case of 
ordinary Newton's law where relativistic effects can't be included directly. 

Using the vector identity, V x (f'A) = /(V x A) — A x (V/), with J = pv, Eq.(6) 
can be written as 

VxJ = Vx(p?7)=p(Vx?7)-?7x (Vp) = , (25) 

which upon using Eq.(4) transforms into 

-l _ v dv , . 

Vxff=- ? x-. (26) 

Thus, Eq.(10) derived from Newton's second law is equivalent to one of the continuity 
equations, viz., Eq(5). Taking the dot product of Eq.(5) with a constant velocity v, we 
get 

v ■ Vp + -irv ■ — = , 
c l at 

*The viscous pressure can be obtained from the viscous force, F = — A?/ 4^ ^> P v = K = —77 V • v 
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which yields 

dp d v 2 
di = di {1 -^ )p - 

According to Lorentz transformation, if the density in the rest frame is p, it will be 
p' = (1 — p in the moving inertial frame. Thus, taking the total derivative, is 
equivalent to taking the partial derivative of the density in a moving inertial frame 
(Lawden, 1968). In terms of the vorticity, the above equation becomes 

v ( dv\ 

* = ~, x -1ST ■ (27) 



c 2 \ dtj 

The vorticity is related to the angular velocity of the fluid (Q) by the relation u = 2 Q. 
This clearly shows that the fluid motion is governed by the continuity equation as well 
as the Newton's equation. 

In a recent paper, we have shown that the magnetic field produced by a moving charged 
particle due to an external electric field is given by (Arbab and Satti, 2009) 

B = -xE. (28) 
cr 

We remark here that there seems to be a resemblance between the vorticity of flow 
and the magnetic field produced by the charged particle. This analogy is evident from 
the fact that B = V x A and uj = V x v. Moreover, Eq.(13) shows that the fluid 
helicity, hf = v ■ u = 0, and Eq.(28) shows that the magnetic helicity, h m — v ■ B — 
(Kikuchi, 2007). Using vector identities, it is obvious from Eq.(ll) that V -u = 0. This 
equation resembles the equation V • B = 0. The former equation implies that vortex 
lines must form closed loops or be terminated at a boundary, and that the strength 
of a vortex line remains constant. Thus, a charged particle creates a magnetic field 
associated with the particle in the same way as the fluid creates a vortex that moves 
with the particle. According to De Broglie hypothesis, a wave nature is associated with 
all moving microparticles. It, thus, seems that vorticity and magnetic fields travel like 
a wave. 

With the same token, Eqs.(27) and (28) suggest that the hydroelectric field is given 

by 

-» dv v -> 

E h = ~l^ u = -xE h . (29) 

This field is generated due to the fluid (mass) motion. Moreover, the Coulomb gauge 
(V ■ A = 0) in electrodynamics is equivalent to the incompressibility of the fluid 
(V • v — 0), however, the Lorentz gauge is equivalent to 

=5 ^ 1 d(v 2 /2) . . 

V- + ? ^P = (30) 

in hydrodynamics^. This latter equation shows that the incompressibility is maintained 
for a steady flow only. When the fluid expands, this merit is lost and the Lorentz gauge 
applies once again. However, from our above definition of hydroelectric field, we must 

*This gauge when multiplied by a constant density p yields the momentum conservation equation. 
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impose the condition that Vx = 0, i.e., x is spatially independent but can depend on 
time, i.e., x = x(t)- The unit of x is J/ kg, or m 2 /s 2 . This defines the energy required 
to move one kilogram of fluid. If we had defined the hydroelectric field in Eq.(29) as 
Eh = — ^ — (where x = "y)> the Euler force in Eq.(10) would become 

F h = -m(E h + vxQ) , (31) 
which is equivalent to Lorentz force in electromagnetism 

F em = q^E e + vxB^j , (32) 

This is a very interesting analogy, since for electron q < so that — m < 0, B Q 
and E e Eh- The complete analogy is tabulated below. We would like to call this 
symmetry an electromagnetohydrodynamics (EMH) analogy. The following table shows 
the analogy between electromagnetic and hydrodynamics. As the electric field of an 
electron points opposite to the force direction, the hydroelectric field (Eh) points oppo- 
site to the direction of flow motion. It is an amazing analogy. Employing this analogy, 
we would like to derive the hydrodynamics laws from the electrodynamics corresponding 
ones. We call this symmetry an electromagnetohydrodynamics (EMH) transitivity. The 
following table shows the analogy between electromagnetic and hydrodynamics. As the 
electric field of an electron points opposite to the force direction, the hydroelectric field 
(Eh) points opposite to the direction of flow motion. 



Theory 


Circulation 


Gauge fields 


Gauge condition 


Helicity 


Electric field 


Electrodynamics 
Hydrodynamics 


B = V x A 
Co = V x v 


if, A 

X = T> * 


V-i+if = 
V-tf+Jrf = 


h e = v ■ B 
fh = v- u 


4 = -f -vx 



It is an amazing analogy. Employing this analogy, we would like to derive the hydro- 
dynamics laws from the electrodynamics corresponding ones. The Maxwell's equations 
are 

VxB = fi J+^-^, (33) 
c z at 

a r> 

V x E e + — = 0, (34) 

V • E e = P -, (35) 



and 



V • B = , (36) 

where p e is the charge density. Employing the EMH analogy, taking the divergence of 
both sides of Eq.(33) and employing Eq.(4) yields 

^ 7? Pm (O r 7\ 
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where ph replaces po and c s replaces c for the EMH, where ph = ^r, where 7 is the 
surface tension and k is the bulk modulus of the fluid under question. Notice that the 
unit of no is Hm" 1 , Eq is FmT 1 and the constant a is kgwT 1 , where ao = t- , which 
may define the resistance (fluidity) of the fluid to flow. Equation (37) is nothing but 
Gauss law in hydrodynamics, where E e — > Eh and the mass density (p m ) replaces the 

2 

charge density (p e ). Besides, Eq.(37) implies the hydrodynamic permittivity, £h = 
Hence, Ph^hC 2 s = 1 (in comparison with PqEqC 2 — !)• Notice however that one can 

2 

relate this constant to Newton's constant by the relation ao = -^q- In this case Eq.(37) 
becomes 

V • E h = AixG Pm , (38) 
which also suggests that G = -p— , which is to be compared with Coulomb constant 



k = 4^-- Therefore, the constant ao is a new fundamental constant. One can therefore 
define the gravitational permittivity, Eh = 1.19 x 10 9 kgm~ 3 s 2 and the gravitational 
permeability as ph = ^r- = 7.25 x 10~ 15 kgm _1 . In such a case the electric field defined 
in Eq.(38) will be gravitoelectric field. This coincides with the value proposed recently 
by Merches and Onuta by considering an ideal gravitomagnetic fluid (Merches and 
Onuta, 1998). The vorticity field intensity of a moving fluid is Hf = — . For the 
Earth u = lCT 14 ^ 1 so that H f = 1.0 kgm~V 2 . Equation (38) can be compared 
with the gravitational field (g) equation due to a static mass distribution of the form, 

V • g — —AirGpm, where p m is a function of space variables only. In this formalism, we 
see that a spinning or orbiting object will generate a hydroelectromagnetic field, where 
the hydroelectric field is equal to the centripetal acceleration. The time dependence of 
the hydroscalar x can be obtained from Eq.(38) using Eqs.(29) and (30) so that 

I|| = ^, x = f (39) 
c 2 at 2 E h 2 

This equation will have direct consequences for cosmological applications where varies 

with cosmic time. For an expanding universe, where a is the scale factor of the universe, 

V • v = 3 (-) = 3H and H is the Hubble constant. Hence, upon using Eqs.(2) and 
(30), Eq.(39) yields the relation p m oc t~ 2 which is a solution of Friedmann equation 
describing the cosmological expansion of a perfect fluid (the universe) (Weinberg, 1972). 
We remark here that Eq.(39) can be applicable to the motion of air inside the lung 
parenchyma. The velocity vector field ( v) can then be obtained by solving Eq.(38) and 
applying the solution in Eq.(29). The solution of Eq.(39) depends on how p m varies 
with time. Using the EMH analogy, Eq.(33) yields 

1 OE 

Vxu = p h J+ ~ 2 -^ • (40) 

This represents Ampere's law in hydrodynamics. Such a result is obtained by (Sulaiman 
and Handoko) in their analogy between electomagnetsism and fluid mechanics [13]. Now 
take the time derivative of Eq.(27) and use Eq.(29) to get 
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Now use Eq.(40) to get 

^ = ?;x(Vxw-/i ft J), J = pv. (406) 

at 

Using the vector identity, V(A- B) = A x (V x B) + B x (V x A) + (A- V)B + (B ■ V) A, 
and Eq.(13), we obtain the equation of vorticity for inviscid flow 

^r + (v- V)u + (u- V)v = 0. (40c) 

It is interesting to note that taking the curl of Eq.(29) and using the fact that 
uj = V x v, one gets 

Vx4+| = 0. (41) 

This is the Faraday analogue of hydrodynamics, where Eh E e and Uh & B. For 
steady flow, Eh = and hence Q = const.. This means that the vorticity is conserved. 
Application of the EMH analogy to Eq.(36) yields 

V • u = , (42) 

which is true using the vector identity V-(Vxu) = V -uj = 0. Equations (38), (40), (41) 
and (42) represent the gravitoelectromagnetic field equations. It is shown by Peng that 
Einstein field equations reduce in the limit of weak field to equations similar to these but 
with negative mass (Peng, 1983, 1990). However, our derivation here doesn't address 
this bizarre situation. This makes our model realistic and physically admissible. Using 
Eq.(40), the fact that u = V x v and the vector identity, Vx(Vxt?) = V(V • v) — V 2 u, 
one obtains 

ig-V^0. (43, 

Equation (43) shows that the velocity vector field, v, is a solution of a wave equation 
traveling at speed of sound. This has the same form as the equation governing the 
electric potential vector A in free space. Hence, instead of solving Euler equation, one 
can solve the above equation to get the velocity. Differentiating of Eq.(43) partially 
with respect to time and using Eq.(29), one obtains 

^#-V^ = 0. (44, 

This equation shows that the hydroelectric field propagates with speed of sound in air, 
and having a source in the density gradient. Hence, the hydroelectric field will arise 
whenever there is any non-uniform spatial distribution of the density of the fluid. It 
is interesting that a static mass distribution can produce a hydroelectric field. It is 
not clear here whether this wave nature of the hydroelectric field is connected with the 
generation of gravitational wave. However, one may assume that any accelerating or 
non-uniform density distribution of the fluid (mass) will induce a hydroelectric field. 
The detection of such a field is to be studied elsewhere. 
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Now, taking the cross product (from left) of both sides of Eq.(40) and using 
Eqs.(41), (42) and (6) the vector identity Vx(Vxu) = V(V-u) - V 2 u), one gets 

This shows that the vorticity vector is a solution of a wave equation traveling at speed 
of sound in air. It is evident that unlike the hydroelectric field the vorticity has no 
source. It is also similar to the evolution of the magnetic field B. Equation (45) can 
also be obtained by taking the curl of Eq.(43). This equation has a spherically symmetric 
solution of the form 

r 

where C(t — r/c s ) is an arbitrary vector of retarded time t — r/c s . The retardation r/c s 
is equal to the time needed for the vorticity wave to pass the distance from the source 
to a given point in space at a distance r. This exhibits the causal behavior associated 
with the wave disturbance. The argument of C shows that an effect observed at the 
point r at time t is caused by the action of the source a distant r away at an earlier or 
retarded time t' = t — r / c s . The time r/c s is the time of propagation of the disturbance 
from the source to the point r. 

According to Maxwell's equations the Poynting vector is given by 

S em = — , (47) 
/Wo 

so that in hydrodynamics, Eq.(47) will become 

S„ = ^, (48) 

fJ>h 

using Eq.(29), Eq.(48) yields 

S h = (p m c 2 s )v. (49) 

This is the same as Eq.(18) with P r = 0. The vanishing of the pressure is due to the 
fact that Vx = 0. Thus, it is remarkable that one can obtain all analogous formulae in 
gravitation by invoking the EMH analogy without deriving them. 



4. Concluding remarks 

We have studied in this paper the consequences of the quaternionic Newton's second law 
of motion. We have derived the energy conservation equation and Euler equation. These 
laws are compatible with the generalized continuity equations. The basic equations 
governing the fluid dynamics are also derived. We have then shown that the energy 
momentum equation is Lorentz invariant. Ordinary Newton's second law is known 
not to be compatible with Lorentz transformation. However, the generalized Newton's 
second laws are compatible with special relativity. Moreover, we have shown that the 
pressure contributes equally to the energy density of the moving fluid as the mass does. 
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Besides, the pressure is the source of the hydroelectric scalar field. Moreover, we have 
found an intimate analogy between electrodynamics and hydrodynamics. With this 
analogy, we derived the magnetohydrodynamics equations from Maxwell's equations. 
We have obtained from Ampere's law, a Gauss's-like law applicable to gravitational 
(or hydrodynamic) system. Moreover, Faraday's-like law and Ampere's -like law in 
hydrodynamics are obtained. In this analogy Euler's force corresponds to Lorentz's 
force. 
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